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One of che oost widely used routines in an algebraic manipulation 
system is a polynomial manipulation package* 1 ' 2,3 *. The crucial 
operation in such routines is the extraction of the Greatest Common 
Divisor (GCD) of two polynomials. This operation is crucial because 
of its freqent use and because it is an expensive operation in regard 
to time and space. 

Experiments by Collins (1) have shown that given two polynomials 
chosen at random, the GCD has a high probability of being 1. Taking 
into account this probability and the cost of obtaining a GCD (some 
GCDs of polynomials of degree 5 in two or three variables can take on 
the order of a minute on the 7094 (,) ), ic appears that a quick method 
of determining whether the GCD is exactly 1 would be profitable. While 
no such complete method is known to exist, a fail-safe procedure has 
been found and is described here. A fall-safe procedure is characterized 
by the fact that when it comes to a decision (in this case that the GCD Is 
1), then the decision is correct. However, the conclusion (i.e. that the 
GCD Is 1) may be true, and the procedure need not arrive at a decision 
regarding it. It is believed that the fail-safe procedure presented here 
(and its extension to the linear case) will arrive at a decision quite 
frequently when the GCD is actually 1. 

We shall first consider the case of polynomials of only one variable. 
The extension to several variables will be made later. 

Let P(x) - a n x" + « n _ 1 x 11-1 + . . . + a Q , * ± integers. n>0 

Q(x) - b n x m + b^j X 1 "" 1 + ... + b Q , b ( integers, m>0 

Let the operation of takiug fan GCD o£ A,B (both for integers and 
polynomials) be represented by (A, 3). 
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Let (P(x), Q(x)) - d(x) and let 

P (x) « d (x) r (x), Q (x) » d(x) 3 <x). 

For any Integer a, Chen 

(P(a), Q (a) ) • d(a) <r<a), 9 (a) > 

Note that (r(a), s<a) ) need not be 1 although 
(r(x), s(x) ) is 1 when the GCD operation is over polynomials. Note further 
that the absolute value of d (a) is always less than or equal to the GCD of 
P(a) and Q(a). 

The plan for the procedure is as follows: We would like to substitute 
a number, a, say, into both P(x) and Q(x) and get the GCD of the results. 
This integer a is so chosen that the leading tera of d(x) will dominate the 
rest of the polynomial, tn fact, we would like the leading term of d(x) to 
be at least twice as large as the rest of d(x). Thus even if the rest of 
the terms of dfx) are of opposite sign froc the leading term, we will be 
able to guarantee that the absolute value of d(a) is at least one-half the 
absolute value of its leading term. If d(x) be of degree k, (k>0), then 
d(a) would be at least Ha k . If <P<a), Q(a) ) is less than ^a k , then d(x) 
cannot be of degree k. If (P(a), Q(a) ) £ ^a , then we do not have a 
decision regarding k. 

In order to find a suitable candidate for a t we will need sone 

information about d(x). Without extensive calculation, it appears that 

we do not know the degree of d(x) and certainly not its coefficients. We 

could , however, find a bound on the zeros of d(x) since such zeros are roots 

of both P(x) and Q(x), 

If we let R be the minimum of R and R , where R and R are bounds for the 

P q P q 

roots of P{x) and Q(x) respectively, then we find the following theorem 

applicable* 

Theorem Let R(>1> bound the absolute value of the roots of d(x), where 

d(x) is a polynomial of degree k>0. Then for x>3kR the absolute value of 

the leading term of d(x) will be greater than twice the absolute value of 

the remainder of d(x). 



3. 



(Note: Thus if we choose a>3kR, where k Is the ainiaua of o and n f 
then we could decide that the GCD Is a constant {i.i, d(x) li of degree 
0) if (P(a), Q(a) ) < H&- «• could decide that the GCD la 1 if we 
previously had taken care of constant CCDs by dividing P(x) and Q(x) 
by the GCD of their coefficients.) 
Proof 

Given that the roots are bounded by R, the condition regarding the 
leading term is hardest to satisfy when the coefficients of the remainder 
of the polynomial are as large as possible and of equal sign. This occurs 
in the polynomial (x+3) . The condition regarding the leading tera can 



now be stated as 



x k > 2 ( (x+R)* - x k ), x>0, R>0 



k>( 



x+R) k 



<£>* X> x + R 



(I) k - 1 



The following leimna will complete the proof. 



It. 



Lonna 



< 3k » k a positive integer 



it >* - i 



Let S 



I L 

i = 2 2 



By the foraula for finite sums of geometric progressions 



,1)- * (1;- -<_(2 r k 
3 m K Z 2 2 



*2' 3 



[f> ' fc - 1 



1 -(f) " fc 



Each terra in S has value <!■ There are exactly k teras in S. Hence S<k 



S = V 3 



-- .-L 



a -(f) " k 



(|)- k 
i - (£) " * 



<3k 



'3* 



(|) * - 1 
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Ve have yet to obtain a complete procedure since we lack & bound 

on the roots of P{x) and Q(xJ. A simple bound for the roots of P{x) 

is n max [ ^i 0$ i£ n - 1] 
an 

This is a crude bound since for (x+RT it yields nB n instead of R. 
Several other easily calculable bounds are known, but we shall leave this 
problem to numerical analysts. 

How let us sum up the procedure defined above. Wc find the minimum 
of the bound on the absolute value of the roots of P(x) and Q(x), Call 
it* R, say. Ve pick an integer, a, say, such that a > 3kR, where k is the 
minimum of the degrees of P(x) and Q(xh If we let g = (P(a), Q(a)) then 
if g< ^ and if ve previously ascertained that the coefficients of P(x) 
and Q(x) had no coercion GCD, then ve can guarantee that d{x) = (P(xl,Q(x})= 1. 
Otherwise no decision is made. 

Extensions to the Basic Procedure * 

Linear GCD 

Let g ■ (?(&)* Q(*)) (where a is chosen by the procedure above). If 

g< x , ve know d(x) ■ 1< Suppose, however, 5 ^ s < 2 at worstr d ^ x ' could 
be linear in X, that is, d(x) ■ ex ♦ d (cj* 0). We propose now to extend the 
scheme defined above in order to determine whether d(x) is linear. If d(x} 
is linear then this extended procedure will find it; if it is not, then d(x) 
must be 1; and then we have extended the range of applicability of the 
procedure for determining whether the GCD is 1. 

As it turns out we can attempt to find all possible linear GCD by 
an extensive search since there are restrictive conditions on c and d- 

Ve know that 

means that e divides f 
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We might be able to avoid testing some or all possible linear 
d{x) by making use of the fact that d(a) divides g. For any pair, 
c,d, such that ac + d does not divide g, ue can conclude that ex + d 
is not a possible d(x). If ac + d docs divide g, then ve shall evaluate 
p C" z) . Q(- 7>- Both of these numbers ttust be zero* If they are, 

d(x) = ex +■ d* If no possible pair of values for c and d works, then 

d(x) = 1. 

Extension to Several Variables 

One can extend the procedure for several variables merely by 
substituting small integers for all but one of the variable and 
determining whether the resulting polynomials have a GCD of 1. If 
they do the original polynomials did alao; if they do not then no 
decision is made. This idea is due to H. >ianove who programed it in 
MATHUB , and existed pri&r to our MSult. It is likely that Hanove's 
idea by itself will yield a very large saving of effort since the 
computation of GCDs grows very badly vhen the number of variables is 
Increased. 

The procedure was attempted on the pair of quartic$ given in (1). 
The evaluation was attempted for values of a between 100 and 150. Most 
CCDs resulted in numbers less than Jga, thus yielding the conclusion of 
a GCD of 1. In the ten cases that did not, the linear case settled the 
matter quickly since ac + d, for possible values of c and d, did not 
divide the GCD with one exception at a - 109* 

For larger pairs of polynomials (say, of degree 30) the evaluation 
will have to be attempted on relatively large numbers (on the order of 
1000), Thus one would have to resort to variable precision arithncric. 
We must consider, however, that the evaluation of the GCD of such 
polynomials by straight forward methods would also require variable 
precision integer arithmetic. 
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It Is hard to estimate, at present, the effect on running time 
that the introduction of a fail-safe aethod will have on the calculation 
of a GCDs. A gain in speed of an order of magnitude or better is not out 
of the question. If the fail-safe technique gains acceptance, then research 
Is likely to yield various methods with differing costs and ranges of 
applicability. 



Evaluation . 

The success of the method depends. In large neasure, on the likelihood 
that the CCD of r(a) and s(a) will be fairly low. Below we give an encour- 
aging result which shows that if r(a) and s(a) were chosen at random, using 
a flat distribution, then the probability that their GCD is greater than 1 
Is less than 0-4. The meaningfulness of this result in practical situations 
is hard to gauge at present. It does, however, give one great hope for the 
success of the method. 

Theorem . Let b be any integer chosen at random from 1,2, ., ,,N 

The probability that an integer less than b is relatively prime to b 
approaches _6 (^0-6) 

as H approaches » 

Proof : The number of Integers less than b and relatively prime to 

It is 

b H (1 - -), where p|b and p Is i prime 

That number can be written as 

b n (1 - - Pr <p|b)), where Pr(p|b) is th* probability 

that p divides b which is 1 if p[b and otherwise. 

Now if we allow b to vary, Pr(p|b) ■ — , that Is, the probability that b 

P i 

is even is \ , that b is divisible by 3 is ^ , etc. 

Thus; the number of integers less than b and relatively prime to It 

becomes , 
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b " (1 " ~ 7> " bn t 1 - 7o>> wh «* p is a prime. 



P P pZ 

But * 
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p* 





n2 " 4 



2 



and 

Thus 

(The theorem used above may be found in any book on the Kieunn Zeta Function) 

This estimate for the probability is already known. The proof above is 
due to W. Henneaan* 

tfe would really be interested in knowing the probability that 
<r(a), s(a)) < e, for small integer values of e. Certain theoretical 
considerations indicate that for e - 10, the probability is greater than 
0.9. However, these considerations are based on certain assumptions re- 
garding the distribution of numbers relatively prime to a given number. 
In experiments we performed these estimates were off by one to two per cent, 
Nonetheless all signs point to the fact that if r(a> and s(a) were chosen 
at random* a high value for their GCD is very unlikely. 
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